Flux flow dissipation in superconductors with short coherence length 
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In superconductors where the coherence length is comparable to the Fermi wave- 
length, the electronic levels within a vortex core are quantized, and separated by en- 
ergies of the order of the superconducting gap. The absence of a continuum modifies 
significantly the energy dissipation due to the motion of vortices. At zero tempera- 
ture, dissipation is suppressed, and the system shows non ohmic behavior. The I-V 
characteristics show two distinct regimes. Below a given electric field, I is strongly de- 
pendent on V. Flux flow resistivity also shows a marked dependence on temperature, 
and reaches conventional values near Tc. 
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Flux flow dissipation is ubiquitous in type-II superconductors. The standard explanation 
of this effect is that, when vortices move, part of the electrical current flows through the core 
(see also 0]). The core region is assumed to be similar to a normal metal. Dissipation 
in the core is described by means of processes indistinguishable from those which take place 
within a normal metal (note also that at fields close to Hc^-, or temperatures close to Tc, 
additional dissipation arises from quasiparticles everywhere in the system 

The basic premise of this argument is that, inside the core, a quasicontinuum of bound 
states exists. The motion of the vortices with respect to the impurities in the lattice gives rise 
to transitions for any velocity or energy. This picture is supported by detailed calculations 
for the low lying states within a vortex core 

Most of the Cu based high-Tc superconductors have coherence lengths at zero temper- 
ature very short, comparable to the lattice constant and to the Fermi wavelength. In this 
regime, usual calculations for the states at the vertex core are not applicable. It has been 
argued that the states within the core should be quantized, and separated by energies com- 
parable to the superconducting gap. Recent observations confirm this picture 0. 

The standard models of flux flow resistivity do not describe microscopically the processes 
which give rise to dissipation, which are taken to be the same as in the normal state. In 
normal metals, there is a continuum of electron-hole excitations. The density of states 
of these modes is proportional to its energy. This assumption suffices to obtain ohmic 
dissipation due to impurity scattering [P,p!0[]. 

In the absence of a continuum, this picture needs to be reformulated. We will assume, 
following the standard work on flux flow resistivity that, when vortices move, some current 
flows through their core. As in conventional superconductors, dissipation takes place through 
the excitation of transitions between states within the core. We study these processes in 
a frame of reference where the vortices are at rest, and lattice impurities are displaced 
towards them. We will consider impurity scattering only, although the scheme can easily 
be generalized to other kinds of disturbances. We assume that impurities are the most 
efficient mechanism at low temperatures. Initially, each impurity is described by the Fourier 
transform of the potential, In the moving frame of reference, the position of a 

given impurity with respect to a static vortex located at the origin is r{t) = fl) + vt, where v 
is the relative velocity of the impurities. Hence, the impurities give rise to a time dependent 
potential, V^e*'^^''e™^*. Dissipation is due to transitions induced by this potential. 

This picture is standard in the study of atom-solid collisions Commonly, the excita- 
tion spectrum of the material, where energy is dissipated, is described in terms of response 
functions. Using this picture for an homogeneous superconductor, we obtain that dissipation 
sets in when fiq^ax^s = where A is the gap of the superconductor, q^ax is the maximum 
wavevector for an excitation at low energies and Vg is the velocity of the superfluid. Setting 
the free electron value Qmax = 2^^, we recover the known result for the critical current in a 
superconductor, in the absence of a magnetic field 0. 

A knowledge of the full response function of the Abrikosov vortex lattice requires lengthy 
and tedious calculations. We will consider first the collision of an impurity with a single 
vortex, and the transitions that ensue. In order to simplify the analysis, we will study a 
two dimensional version of the problem. The vortex defines the z axis. For very anisotropic 
superconductors, hke many high-Tc materials, a 2D picture is fully appropiate. In general, 
core states give rise to bands along the z direction. Transitions are characterized by the 
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wavevectors, k^, k'z ~ kp, of the states involved. If the impurity potential is sufficiently 
short range, we can ignore the dependence of the transition rates on these momenta. If we 
neglect as well the dispersion of these one dimensional subbands, and describe each of them 
by a typical energy, e^, the problem becomes again two dimensional. The analysis to be 
discussed below applies, with the only inclusion of a prefactor kpaz, where is the range of 
the impurity potential along the z direction. In any case, the most dramatic deviations from 
standard flux flow dissipation will be shown to be originated by the existence of a sizeable 
energy gap within the vortex core. This is a general feature of superconductors with a short 
coherence length. 

A simplifled sketch of the impurity vortex collision is shown in the flgure. Within second 
order perturbation theory, the matrix element of the impurity potential and the states 
involved in the transition are required. The quasiparticle states induced at the vortex core 
have a smooth envelope, with a characteristic localization length of the order of the coherence 
length of the superconductor, ^, and rapidly varying oscillations, characterized by kp^, where 
kp is the Fermi wavevector []6|,[7|,p!2[| . Thus, the spatial dependence of these wavefunctions, 
within a plane normal to the vortex, is like ~ ^-ig-H/Cg^'^Fr^ addition, shows an 
angular dependence. Its influence in dissipative processes, however, is much weaker than 
the modulations discussed above. 

For convenience, we will assume that the range of the (2D) impurity potential is shorter 
than the length scales mentioned above. Then, we can neglect its q dependence, and write it 
as J2g Va^e'^'^^, where V is the strength of the potential, and a is a length of the order of the 
range of the impurity potential. We focus on a single impurity, which follows a trajectory 
r = ro + vt. The vortex is located at the origin of coordinates. The transition amplitude 
between two vortex states with an energy splitting e^j is given, within second order 
perturbation theory, by: 



In terms of the transition amplitude A, the energy dissipation per unit time and per unit 
vortex length can be written as: 



where the summation is over all intravortex transitions, and an integral over all impact 
parameters, ro, is required. This integral is cut off by the smallest separation, Vmin, possible 
between the impurity and the vortex center. This distance is of the order of the lattice 
spacing, which we can take similar to kp^. n is the (3D) impurity concentration. 

The integral over time in eq. (1) is determined by times such that t < tmin^/^, where we 
are assuming that tmin is the shortest time in the integrand. Its behavior, as function of ro 
shows two different regimes, depending on the value of hv/{ero). If hv/e <^ tq, the integral 
has an exponential dependence on e/{hvmm{kp^,$,,ro). Due to the existence of a short 
distance cutoff in the integral over ro, this can be the relevant regime for superconductors 
with shoet coherence length, as discussed below. In that case, the energy dissipation can be 
written as: 
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It is worth noting that, on the other hand, equation (2) gives rise to the standard ohmic 
dissipation for a metalhc density of excitations. For an extended, normal vortex core, we 
can replace the sum over possible transitions by an integral. The density of low energy 
excitations in a normal metal goes as A^(e) ~ ^/Ep. Inserting this expression in equation 
(2), we find that dE/dt ~ f^, which leads to ohmic behavior. 

Equations (2) and (3) give the energy dissipation due to a single vortex. The dissipation 
per unit volume can be obtained by multiplying by the (two dimensional) density of vortices, 
~ H / {H c^^"^) 1 where H is the applied field. 

In order to relate the dissipation to the I-V characteristics of the material, we will assume 
that the voltage drop is due to the vortex motion ||13|. Then, the electric field within the 
sample is \E\ = v\H\/c. The dissipation per unit volume is given by aff\E\'^, where c// is 
the flux flow conductivity. 

To our knowledge, there are no direct measurements of the vortex velocity in a type- 
II superconductor. For magnetic and electric fields typical for high-Tc materials, H ~ IT, 
E ~ IV/cm, we find that f/c ~ 10~^ — 10~^. This result implies that, when the lowest energy 
splitting within the vortex, e ~ A/lkpC,)' is e ~ hv/r 

min ~ 2meV, the effects considered here 
become relevant (we are taking r^m ~ lA). This value is below the estimated gap of high-Tc 
superconductors, A ~ 3 — Sfc^Tc ( assuming conventional s-wave pairing). This estimate 
can be reduced for larger magnetic fields or smaller voltages. 

The dependence of the transition rates on v, for sufficiently low values of v, is determined 
by the exponential in eq. (3). Expressing the dissipation as jE, and using the fact that 
E oc V, we conclude that I oc exp —{Vo/V), where Vq is a threshold voltage. Its value is such 
that the local electric field is given by Eq ~ eVminB/ (he). For the parameters mentioned 
above, this field is 1 V/cm. For fields lower than this value, dissipative processes are strongly 
suppressed. Note that this effect should be more easily observable at high magnetic fields, 
as the velocities of the vortices are smaller, and the motion is less effective in inducing 
transitions. This decrease in the vortex velocities more than compensates the increase in 
the number of vortices induced by the magnetic field. Thus, flux flow dissipation, for a 
fixed voltage, should decrease with increasing magnetic field, in sharp contrast with the 
Bardeen-Stephen model [Q. 

The anomalous dissipation in superconductors with short coherence length can be masked 
by pinning effects. Pinning energies also show unconventional dependence on applied cur- 
rents, magnetic fields and temperature in high-T^ materials. The effects reported here in- 
fluence the flux dynamics mostly at low temperatures (short coherence lengths), which also 
enhance the pinning of vortices. High magnetic fields, on the other hand, reduce pinning, 
and require lower vortex velocities to generate a given voltage drop. Thus, we expect that the 
suppression of vortex viscosity will be most effective in modifying the properties of high-Tc 
superconductors at low temperatures and high magnetic fields. 

In many situations, flux creep determines the magnetic properties. The dynamics of 
vortices are described in terms of an exponential dependence on an activation energy, and 
a characteristic frequency Jl^, which describes the motion of the vortex within each po- 



tential well. This frequency, when the motion of the vortices is overdamped, depends on 
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the viscosity. So far, its dependence on magnetic field and temperature has been neglected. 
The analysis reported here suggests that this prefactor can also have a strong dependence 
on fields and temperatures. As function of temperature, we expect a crossover from a high 
temperature regime, where the motion is overdamped, and this attempt frequency is deter- 
mined by the viscosity alone, to a low temperature regime, where the viscosity is negligible, 
and the attempt frequency depends on the mass of the vortex. 

Our results may help to explain anomalous effects in the ac response of high-Tc supercon- 
ductors in megnetic fields |T^. The strong dependence of / on V, for low voltages, described 
above may help to explain a variety of experiments, which have been fitted to J ~ V"" laws, 
with a ^ 1 [|T6|. The absence of dissipative mechanisms, at low temperatures and high 
magnetic fields may also suppress the magnetic relaxation rate [jl^. Magnetic relaxation 



involves the exchange of energy from the magnetic field into thermal energy. Our estimates 
for the relevant vortex velocities are also very close to those of interest in the vortex tunnel- 
ing regime . A change in the dissipation mechanisms can induce drastic modifications in 
the tunneling rates, as they depend exponentially on the viscosity. 

It is interesting to speculate with the behavior of a type II superconductor in the limit 
of very large energy splittings within the vortex cores. Then, conventional dissipative mech- 
anisms are frozen out, and, in the absence of pinning centers, vortices will move without 
friction. The motion of vortices gives rise to time dependent magnetic fiuxes and induced 
electric fields, as discussed above. Thus, a voltage drop accross the superconductor will be 
generated. On the other hand, the absence of dissipation impedes the existence of currents. 
Unlike a superconductor withouth an applied magnetic field, such a system should be able 
to sustain a voltage drop, without induced currents. This unusual behavior, which closely 
resembles that of an insulator, will be masked by the existence of pinning, which can localize 
the vortices. 

The analysis presented so far relies on the applicability of the BCS theory to supercon- 
ductors where ^ ~ kp^. The limit of a strongly coupled superconductor, however, loosely 
resembles an array of weakly coupled superconductiong islands. In these systems, local su- 
perconductivity can appear at temperatures well above those where global coherence sets 
in. This is similar to the expected behavior of a strongly coupled superconductor, where 
preexisting pairs are supposed to exist above T^. Dissipation at low temperatures in arrays of 
Josephson junctions is being studied extensively, both experimentally [0 and theoretically 
p!9| . It seems confirmed that, in the absence of nonsuperconducting, resistive channels, dis- 
sipation is strongly reduced at low temperatures. This result can be understood on general 
grounds, since the superconducting system has no low energy excitations below the plasma 
frequency. 

In conclusion, we have extended the standard theory of fiux fiow dissipation to super- 
conductors where the level splittings, e, between the states within the vortex core cannot 
be neglected (e ~ fiv/^, where v is the velocity of the vortices). In this regime, dissipation 
is strongly suppressed. The response of the superconductor is non-ohmic, and I shows a 
marked dependence on V. The resistence is anomalously large. At a given voltage, the 
system undergoes a crossover into a more conventional flux flow regime. 

Simple estimates show that this analysis may be relevant to high-Tc superconductors, at 
sufficiently low temperatures, so that kp^iT) ~ 1. The currents should be large enough, in 
order to avoid pinning effects. It would be desirable to explore this regime of temperatures, 
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magnetic fields and currents |2T|, in order to check the results outlined here. 

This work has been partially supported by CICyT (grant MAT91-0905). We acknowledge 
illuminating discussions on collision induced transitions and dissipation with Prof. P. M. 
Echenique. 
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FIGURES 

FIG. 1. Schematic picture of a vortex-impurity scattering process. See text for details. 
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